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Nonlinear phase shift from photon-photon scattering in vacuum 
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We show that QED nonlinear effects imply a phase correction to the linear evolution of electro- 
magnetic waves in vacuum. We provide explicit solutions of the modified Maxwell's equations for 
the propagation of a superposition of two plane waves, and calculate analytically and numerically 
the corresponding phase shift. This provides a new framework for the search of all-optical signatures 
of photon-photon scattering in vacuum. In particular, we propose an experiment for measuring the 
phase shift in projected high-power laser facilities. 
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Introduction.- Quantum Electrodynamics (QED) pre- 
dicts that at photon energies well below the elec- 
tron rest energy, photon-photon collisions can still be 
produced through the interchange of virtual electron- 
positron pairs P, This nonlinear interaction modi- 
fies Maxwell's equations for the average values of the 
electromagnetic quantum fields Q and affects the prop- 
erties of the QED vacuumQ. During many years, the 
search of these effects has been restricted to projected 
particle physics experiments with accelerators. However, 
photon-photon scattering processes in vacuum will be- 
come testable at energy densities achievable with ultra- 
high power lasers in the near future 

The realization of these ultra intense photon sources 
began with the discovery of chirped pulse amplification 
(CPA)Q in the late 80's and optical parametric chirped 
pulse amplification (OPCPA) Q in the 90's. These tech- 
niques opened the door to a field of research in the 
boundary between optics and experimental high-energy 
physics, where lots of novelties are expected to come in 
the next years. In fact, several recent works propose dif- 
ferent configurations that can be used to test the nonlin- 
ear optical response of the vacuum, e.g. using harmonic 
generation in an inhomogeneous magnetic field [8[, QED 
four- wave mixing^, resonant interactions in microwave 
cavities [10}, or QED vacuum birefringence which can 
be probed by x-ray pulses among othersjl^. 

In the present Letter, we show that QED nonlinear ef- 
fects in vacuum imply a phase correction to the linear 
evolution of crossing electromagnetic waves. We provide 
explicit numerical and analytical approximate solutions 
for the propagation of a superposition of two plane waves. 
This result allows us to calculate the corresponding phase 
shift, providing a new framework for the quest of signa- 
tures of photon-photon scattering. In particular, we sug- 
gest an experiment for measuring the effect of nonlinear 
vacuum in projected high-power laser facilities like the 
European Extreme Light Infrastructure (ELI) project 14 
for near IR radiation. 

Model and equations-. Let us begin by writing the for- 
mulae that have to be used instead of the classical linear 
Maxwell equations, including the terms which come from 



QED effects in vacuum. The corresponding Lagrangian 
density in terms of the electric and magnetic fields E and 
B, was derived in the 30's by Euler and HeisenbergQ: 



C = C + £C Q = Co + £ 
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being 



C n = f£ ( E z - c z B' 
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the linear Lagrangian density and eo and c the dielectric 
constant and the speed of light in vacuum, respectively. 
As it can be appreciated in Eq. ([1]) , QED corrections are 
introduced by the parameter 
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This quantity has dimensions of the inverse of an energy 
density. This means that significant changes with respect 
to linear propagation can be expected for values around 
|££o| ~ 1) corresponding to beam fluxes with electromag- 
netic energy densities given by the time-time component 
of the energy-momentum tensor 
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d t A - C > 2/£ ~ 3 x 10 29 J/m 3 



(4) 



While such intensities may have an astrophysical or cos- 
mological importance, they are not achievable in the lab- 
oratory. The best high-power lasers that are being pro- 
jected for the next few decades will be several orders of 
magnitude weaker, eventually reaching energy density of 
the order p ~ 10 23 J/m 3 [5]. Therefore, we will study here 
the "perturbative" regime, in which the non-linear cor- 
rection is very small, "C 1. As we shall see, even 
in this case measurable effects can be accumulated in the 
phase of beams of wavelength A traveling over a distance 
of the order A|££o| -1 - Thus, current sensitive techniques 
could be used to detect traces of QED vacuum nonlin- 
earities. 
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Once the electromagnetic fields are expressed in terms 
of the four-component gauge field — (A , A) as B = 
V A A and E = -cVA° - ^ 
are given by the Variational Principle 



cVA° — the equations of motion 



ST 



o, 



(5) 



where T = J Cd x is the QED effective action. Instead of 
studying the resulting equations for the fields E and B, 
that can be found in the literature [H, [HI, for the present 
purposes it is more convenient to consider the equations 
for the gauge held components A^ 1 . In general, these four 
equations cannot be disentangled. However, after some 
straightforward algebra it can be seen that they admit 
solutions in the form of linearly polarized waves, e.g. in 
the x direction, with A = and A = [A, 0, 0), provided 
that: i) the held A does not depend on the variable x 
(a transversality condition) and ii) A(t, y, z) satisfies the 
single equation: 



d^A + £e (d^Ad v Ad v A + 2d^Ad»Ad v A) = 0, (6) 

where we have used the convention = 
diag(l, — 1, — 1, — 1) for the metric tensor. Here- 
after, we will restrict our discussion to this case of 
linearly-polarized solution. The orthogonality relation 
E • B = are then automatically satisfied, and the ef- 
fective Lagrangian Eq. (P) reduces to C = Co (1 + ££o)- 
Note that the plane- wave solutions of the linear Maxwell 
equations, such as .4 cos (k • r — uit), where A is a 
constant, k = (0, k y , k z ) and lo = c|k|, are still solutions 
of Eq. ©. However, we expect that the non-linear 
terms proportional to £, due to the QED correction, will 
spoil the superposition principle. 

Variational solutions-. Let us consider two plane 
waves, for simplicity having the same phase at the space- 
time origin, having wave vectors (0, q, k) and (0, —q, k) 
respectively, and angular frequency lu = \Jk? + q 2 . Any 
of them, when taken alone, would be a solution of both 
the linear and non-linear equations. However, the super- 
position 



A(t,y,z) 



A 



[cos(fcz — Lot + qy) + cos(fcz — Lut — qy)] 



Acos(qy) cos(fcz — u>t), 



where A is a constant, would only solve the linear equa- 
tions of motion. In our perturbative regime, we can ex- 
pect that the small non-linear correction will progres- 
sively modify the form of A(t, y, z) as the wave proceeds 
along the z direction. We will therefore make the ansatz 

A = Acos(qy) [ a(z) cos(fcz — cut) + (3(z) sin(fcz — ut) ] , 

allowing for the generation of the other linearly- 
independent function sin(fcz — tot) (we will take a(0) = 1 



and (3(0) = 0). On the other hand, we neglect the possi- 
ble generation of a reflected wave depending on kz + ut, 
which can be expected to be a smaller correction in this 
perturbative regime. As we will discuss below, in our per- 
turbative regime the effects of the possible y-dependencc 
of a and (3 are negligible for the traveling distances in 
the z direction that we will consider. 

According to the Variational Method, we require that 
the functions a(z) and (3(z) correspond to a local mini- 
mum of the effective action T, after averaging out y and 
t as follows: 



whose minimum corresponds to the equations ST /5a = 
y ST/S/3 = 0. After a straightforward computation, these 
two equations can be written as 



^- + k0 + aq A a{(3 2 + « 2 ) - ^ [2 (9k 2 - Aq 2 ) (3 2 (9) 

+54k(3a' + 27a' 2 + 6k 2 a 2 - T8kaj3' + 90 2 ] a" 

+^ \9kaa' - 9a' - 9k(30 + 2 (3k 2 - 2q 2 ) a0] (3" 
8 

+ ^ [_9 fc/ 3'3 _ 3fc 2 ac/ 2 + ( 15fc 2 _ g? 2) a/? /2j 



-- \9ka' 2 - 16kq 2 a 2 - \%kq 2 + 2 (9fc 2 - 4 9 2 ) a'0] 
8 



= 0, 
and 



Y~ ka ' + «9 4 /?(« 2 + I3 2 ) - ^ [2 (9k 2 - 4q 2 ) a 2 (10) 

-54ka0 + 27/3' 2 + 6k 2 + I8k(3a' + 9a' 2 ] (3" 
+ | [-9kf3f3' - 90a' + 9kaa' + 2 (3k 2 - 2q 2 ) (3a] a" 

ka' 3 - 3k 2 (30 2 + (15k 2 - 8q 2 ) (3a' 2 ] 

-9k0 2 + I6kq 2 (3 2 + I6kq 2 a 2 + 2 (9fc 2 - 4? 2 ) 0a] a' 



a 

"8 
= 0. 



where a = £eoA 2 /2. It is now convenient to express the 
parameter a in terms of the energy density 

P = Too (11) 
= ^(E 2 + c 2 B 2 ) + 1 2 (E 2 - c 2 B 2 ) (3E 2 + c 2 B 2 ) . 



n 2 

In our perturbative regime, this gives a ~ jjz~£,P with a 
very good approximation. 

Numerical simulations and approximate analytical 
solution-. The result of the numerical integration of Eqs. 
© and (I10p is shown in Fig. 1 for a choice of parame- 
ters that may be accessible at future facilities[5|, namely: 
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p = 4 x 10 23 J/m 3 , A = 5 x 10~ 7 m, k = q = = 

0.89 x 10 7 m -1 , giving a = 3.4 x 10~ 20 m 2 . The two en- 
veloping functions show a sinusoidal behavior, with an 
oscillation length ~ 0.26m. In particular, we find that a 
100% change in phase is accumulated after the distance 
Az ~ 13cm. 




z(m) 

FIG. 1: Numerical solution of Eqs. © and iflu) for fe = g = 
0.89 x 10 7 m- 1 and a = 3.4 x 10~ 20 m 2 . The continuous curve 
represents the function a(z), the dashed curve the function 

This numerical result also suggests the viability of an 
analytical approximation. In fact, the numerical oscil- 
lations show that each order in derivation of a and 
corresponds to a suppression factor ~ 27r/0.26m _1 , 6 or- 
ders of magnitude smaller than k and q. We can then 
neglect the second derivatives in Eqs. ([9]), (fTOf . More- 
over, in our perturbative regime a can also be considered 
an expansion parameter, therefore we can also neglect the 
first derivatives when they appear multiplied by a. As a 
result, Eqs. |9j), (fTTJ|) can be approximated as follows: 

k0 + aq 4 a(a 2 + (3 2 ) = 0, (12) 



- ka' + aq i (3{a 2 + /3 2 ) = 0, (13) 
whose analytical solution is 

a(z) — cos (x z ) ; an d f3{z) = — sin (xz) , (14) 
where 



X = aq^k = ie a A 2 q 4 /2k ~ 2^-%p. (15) 



of the parameters as in Fig. 1. This is an a posteriori 
justification for the analytical approximation. 

By substituting in Eq. ([7]) and after some trigonome- 
try, we find that our approximated analytical solution of 
the non-linear equations for the electromagnetic field in 
the vacuum can be written as 

A(t, y, z) — Acos(qy) cos(kz + x z ~ <*>i). (16) 

Therefore, the effect of the non-linearity is to change 
the phase of the wave, with respect to the linear solution, 
by a term that increases linearly with the distance. This 
can also be interpreted as a change in the z component of 
the wave vector, which becomes k z = k + x, so that the 
dispersion relation is modified to u> — cy/ (k z — x) 2 + q 2 ■ 

To conclude this section, let us discuss now the validity 
of the approximations that we have made. The Varia- 
tional Method that we have used is expected to provide 
very good results whenever the class of test functions, 
given by Eq. (0 in our case, is a reasonably good choice. 
Since the latter was motivated by perturbative considera- 
tions, the whole approximation will be justified whenever 

2 4 

the correction x <C k, i.e. whenever <C 1. This 

is guaranteed by the dispersion relation and the fact that 
£p -c 1 in our regime. It is easy to see that in this case 
the correction to the energy density dependence on k and 
q is also very small. 

The other approximation that we have made was ne- 
glecting the y-dependence of the enveloping functions a 
and (3. For an a posteriori test of this hypothesis, we have 
corrected our solution Eq. (fT75)) allowing for the further 
envelop functions j(y, z) and a(y, z). We have then used 
as a new ansatz for the Variational Method the potential 
A = Acos(qy)('j cos(k z z — tot) + asm(k z z — u>t)), where 
k z — k + x = k + aq^/k from Eq. (fTS"]) . 

We have performed a first average over the fast vari- 
ation in y in the action due to the trigonometric depen- 
dence on the product qy, and we have minimized the 
action. After keeping only the first order in the ex- 
pansion parameter a, and neglecting the terms involving 
the second derivatives (which imply a slower variation 
as discussed above), we have got the equations d z a + 
XI (l 2 + cr 2 - l) =0, and d z j - x° {l 2 + o 2 - l) =0. 
These equations, which do not involve d y j and d y a at 
the first order, are solved by j(y, z) = 1 and a{y, z) = 0, 
therefore we conclude that at this order our previous so- 
lution, Eq. (fT6|) . is not modified. This justifies the ap- 
proximation of neglecting the possible y-dependence in 
the first instance. 

Comparison with the optical Kerr effect-. Since it is 
theoretically known that the vacuum shows birefringence 
as in the DC Kerr effect [4], it is interesting to compare 
our result with the optical (AC) Kerr effect that also 
arises in matter. In fact, let us consider a Kerr medium, 
characterized by an effective nonlinear refractive index of 
the form: 



This result fairly coincides with our numerical solution, 
giving x = 24m" 1 and 27r/x = 0.26m for the same choice 



n = n + n 2 J, 



(17) 
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where no is the linear refractive index, n 2 the Kerr coeffi- 
cient and I the irradiance of the beam. In such a medium, 
for propagation through a distance dz, the wavefront 
phase will be modified by an amount: 



d$ = —n 2 Idz. (18) 
A 

Since p = I/c, Eqs. (fT8|) and (TT5]) show that our con- 
figuration of two crossing waves in vacuum undergoes the 
same phase shift as a single plane wave in a Kerr medium 
having uq — 1 and a Kerr coefficient n 2 = 2c£q 4 /(uj 2 k 2 ). 
Choosing e.g. q — k — w/(v2c), this Kerr coefficient 
would be n 2 = £/(\/2c) w 10~ 38 m 2 W _1 . However, in 
spite of this analogy, it is important to note that such an 
n 2 cannot be interpreted as a Kerr coefficient for the vac- 
uum. In fact, strictly speaking the vacuum does not show 
the usual AC Kerr effect. In a Kerr medium, the phase 
shift (|18p is found even for a single plane wave propa- 
gating along the z-direction. In vacuum, such a single 
plane wave, corresponding to q — in Eq. (|15[) . would 
propagate with x — 0, i.e. without any phase change, 
exactly as in the linear case. Ultimately, this is due the 
fact that the cross section for photon-photon scattering 
vanishes for parallel momenta. 

Proposal of an experiment-. Our previous discussion 
suggests that we might test the non-linear properties of 
the QED vacuum by using sensitive experiments in which 
the key point is the ability of measuring small phase 
changes in a laser beam. In our configuration, a laser 
pulse is divided in three beams of the same intensity and 
two of the resulting beams are focused independently. 
The trajectories of both rays cross at the focus point 
with an angle 6. As a result, the central part of each dis- 
tribution has acquired a phase shift A<I>. In a experiment 
corresponding to the parameters of the ELI project in its 



first step we have pulses of wavelength A = 800nm, in- 
tensity / = 10 29 Wm~ 2 and duration r = 10/s which are 
focused in a spot of diameter d ~ 10/im. Either from Eq. 
(|15|). or equivalently from Eq. (TT5)) . this results in a phase 
shift A<E> « 10 -7 rad, which can be resolved comparing 
with the third beam which was not exposed to the ef- 
fects of QED vacuum. Current techniques like spectrally 
resolved two-beam coupling, which can be applied for 
ultrashort pulses [Ttij. can be used to this purpose. Com- 
paring to other alternatives like x-ray probing of QED 
birefringence, our system does not need an extra free 
electron laser and the power requirements of the system 
are only one order of magnitude higher. Moreover, the 
measurement of the ellipticity and the polarization rota- 
tion angle in birefringence experiments is not yet possible 
with current technology. Other techniques like four-wave 
mixing processes @ require the crossing of at least three 
beams, with the corresponding alignment problems and 
the rest of the requirements are similar to our proposal. 

Conclusions.- In conclusion, we have calculated the 
phase shift arising from propagation of ultra-intense ra- 
diation in vacuum, and shown that it could be measured 
in the first step of ELI facility under construction. We 
consider that the present work could serve as a starting 
point for the quest of other nonlinear optical phenomena 
that may arise in ultra-high power laser beams propagat- 
ing in vacuum. 
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